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Abstract 

We calculate the one- loop QCD correction for inclusive jet cross section and Drell-Yan 
process in a general low-energy effective Lagrangian for composite quarks and leptons. 



1 Introduction 

Recently the CDF group at Fermilab reported a significant excess of one-jet inclusive production 
at high pt over the standard QCD prediction [|[]. The observed inclusive jet cross section is 
in excellent agreement with theory, but begins to deviate from the QCD prediction around 
Pt = 200 GeV, and its central value becomes as large as twice the theoretical prediction at pt 
above 400 GeV. If this discrepancy between theory and experiment survives more stringent tests, 
and arises not from the uncertainties in the QCD parameters such as the parton distribution 
functions but from a genuine new physics, one possible new physics explanation would be 
existence of a quark substructure. 

A substructure in quarks gives rise to four-fermion contact interactions at small energies 
compared to the compositeness scale A via constituent exchanges, and this induces a correction 
of order s/A 2 to the QCD prediction of jet production f2|. The correction is negligible at 
small energies, but becomes significant at high p?. This behavior agrees qualitatively with the 
observed inclusive jet cross section. We assume here that only quarks are composite and gauge 
fields are elementary. 

The CDF fit of the data using the tree level amplitudes from the effective Lagrangian 
by Eichten, Hinchliffe, Lane and Quigg (EHLQ) M with SU{2) L doublet quarks gives the 
compositeness scale A w 1.6 TeV. To go beyond the tree level analysis of the data, we need 
the QCD one-loop corrected amplitudes for q^qj — > q^qj. The leading QCD correction to the 
amplitudes arising from the QCD interaction has been known || 0, and so only the one-loop 
correction to the terms arising from the four-fermion contact interactions needs to be computed. 

In this paper, we calculate the one-loop QCD correction to q^qj — > q^qj in the EHLQ effective 
Lagrangian, using the framework of Kunszt,et.al. for one-jet inclusive cross section || [/J, and 
also discuss QCD corrections in Drell-Yan process. Section 2 through 6 is devoted to QCD 
corrections for the inclusive one-jet cross section. In sec. 2 and 3, we review the EHLQ 
Lagrangian, and give the squared amplitudes at tree-level, and in sec. 4 discuss the ultraviolet 
divergence and summarize the short distant effects of loop corrections. In sec. 5 we briefly review 
the method by Kunszt,et.al. for one-loop inclusive jet cross section, and give our result in sec. 6. 
Details of QCD calculation may be found in the appendix. Finally in sec. 7 we discuss QCD 
corrections in Drell-Yan process. 



A typical term of the helicity conserving effective interactions of composite quarks at low 
energies compared to the compositeness scale can be written in the form of current product 
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where rjfai, A) and are renormalized effective coupling and generic quark current respectively. 
The constituent exchange between currents is represented by D M ^(x,A) which is assumed to 
satisfy 

ZV(x,A) = A 2 v ^(iA), (2) 
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x,A) - g^S^(x) for xA > 1. (3) 



The A in D^ u is a cutoff that determines the interaction range of the constituent exchanges. The 
relation between rj(fi, A), and the corresponding bare quantities depends on the fundamental 
dynamics at the compositeness scale. However, this model dependence does not cause any 
problem in calculating QCD corrections at low energies because any ambiguity arising from the 
lack of knowledge on how the currents and couplings are renormalized can be absorbed in the 
coupling r] which is supposed to be determined experimentally. With (|), ([!]) becomes at tree 
level 

f?^- (4) 

The most general helicity conserving, SU(3) C x SU(2)l x symmetric low-energy ef- 

fective Lagrangian of up and down quarks by EHLQ is 



( q 2 \ ( T a T a 

£ehlq = 9 2 l^^i^J Y] q L ^q L q L YqL + ViQLl^Y qL ^ L ^Y qL 
+VuqLl fl qLU R 'J tl UR + ri d q L ^q L d R Yd R 

+m u q L ^—q L u R n ti —u R + rimLi^—qLdRi 11 —d R 

+VuuU R 'j f *u R u R 'y> J 'u R + ri dd d R ^d R d R ^d R 

+r]udU R YuRd R Yd R + ij' ud UBr/' i d R d R 'y ti URj , (5) 

where qi = (uL,d L ). We inserted in (|5|) the strong coupling g 2 explicitly to make the tree 
amplitudes of QCD and contact terms be formally in the same order in the QCD coupling. For 
convenience, in the rest of the paper we absorb the factor 

g2 ° (6) 



2g 2 A 2 



into the coupling 77's. We also assume here all quarks are massless . Then because of the 
SUl{2) symmetry, there are only seven independent helicity amplitudes for qiqj — ► q^qj. They 
are: u L d L -> u L d L ,u L u L — ► u L u L ,u R d L -> u R d L ,u L d R -> u L d R ,u R u R -> u R u R ,d R d R -> d R d R , 
and u R d R — > u R d R . The amplitude for d L d R — ► d L d R , for example, is identical to that of 
Uhd R — > uid R because of the SU{2)l symmetry. In the following, we calculate these seven 
amplitudes to one-loop order in QCD. 
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3 Tree amplitudes 



The tree level amplitudes for the helicity channels in the EHLQ effective action are given in the 
appendix. Note that we follow the notation in ref. 0, |8j for the helicity amplitude and spinor 
algebra. 

The squared amplitudes — color and spin averaged — for quark channels are 

2 



\A(ud — > ud)\' 
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A{ud — > ud) 
i 



^ S ^ U + u 2 (vl + Vd + ^ivL + vld) 



2 11 2 

+s 2 (Aril + -tjoVi + j^vl + ^VudV'ud + vld + Vud) 

8 s 2 Au 2 

+ gy(^l + Vud) + gj(*« + Vsd) 



(7) 



\A{uu — > uu) 



\A{uu — > uu)Y 
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4 (s 2 + u 2 | s 2 + t 2 2s 2 ^ 
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t 2 



3tu J 



4 (s 2 s 2 " 



+ ^-(16% 2 + 8 W i + vl + ^Vlu) + 2(n 2 + t 2 )( V l + ^ 2 J 



and using the crossing symmetry 



A(ud — > w<f) = |^4(tZc£ — > mc/)| = |^4(-u(i — > uc 



u 



A(uu — > dd) = A{dd — > mm) = |.A(tt<i — > ud)| 2 (s — > — ► s, u — > t), 
|^4.(wu — ► uu)| 2 = — > im)| 2 (s <-» m), 



with 



(8) 

(9) 
(10) 
(11) 

s = (pi + p 2 ) 2 , t = (pi - P4) 2 , w = (pi - P3) 2 - (12) 
For the channels involving d-quarks only, we can obtain the squared amplitudes by replacing 
Vuu,Vsu with r)dd,V&d respectively in (|),(0). 
Putting back the factor 

9ll±[_ , s 

2a s A 2 [ } 

into each 77's in the above equations and keeping only r] Q we can recover the formulas in O] with 

g 2 /An = 1. (14) 

Note that we have corrections for typos in (8.13), (8.15) in B. 



4 Ultraviolet divergence 



The one-loop Feynman diagrams for qiqj — > in the EHLQ effective Lagrangian are given in 
Fig.l. In the massless limit, the one-loop self energy diagrams for fermions in the dimensional 



regularization vanish |TIJ, and (g),(h) in Fig.l are the UV counter terms for current renormal- 
ization arising from the fermion self-energy diagrams and (a),(b). (g),(h) are required only for 
color octet currents, and the counter terms for the conserved, color singlet currents vanish. We 
do not include in our calculation the one- loop diagram in Fig. 2(b) because its finite part coming 
from the small momentum region (~ y/s) is suppressed by a factor s/A 2 relative to those in 
Fig.l. This can be easily seen from the fact that the four-fermion-gluon vertex in the diagram 
is given by the effective Lagrangian 

~ -j^fabc Jl JlF^ v (15) 

represented in Fig. 2 (a). Here J® denote color octet currents. The contribution from the large 
momentum region (3> y/s), which is dependent on the fundamental dynamics at the compos- 
iteness scale, is independent of external momenta, and renormalizes only the couplings of the 
contact terms. Since we are not interested in how the 77's are renormalized, we can completely 
exclude this diagram. 

There are also penguin diagrams (Fig.3). Although they are not one-loop QCD corrections, 
it is easy to see that they induce form factors in the quark-gluon vertices that are formally of 
same order of magnitude as the QCD one-loop corrections. The form factors induced by the 
penguin diagrams assume the form 

where C (ln(— q 2 )) is a model dependent function of 0(1). In the last step we replaced the 
function C with its average value C in the momentum range of interest. Thus penguin diagrams 
introduce new free parameters in the amplitudes. This form factor effect from penguin diagrams 
may be combined with that in the vector boson propagators 0, [| . 

Now in general, (a)-(f) have ultraviolet divergences as well as soft and collinear divergences. 
For one-jet inclusive cross section, the soft and collinear divergences are cancelled by those from 
2 — > 3 process, which will be reviewed in more detail in the next section. The UV divergences 
in (a), (b), as mentioned before, are cancelled by the counter terms (g),(h). We assume the 
counter terms are given in the MS scheme. The scheme dependence of the counter terms 
is absorbed in the coupling 77's to render the physical amplitudes scheme independent. The 
UV divergences in (c)-(f) arise from the approximation in (Q). If we insert D llv {x,A) between 
the currents, the diagrams would be finite with logarithmic terms of order a s log(A) from the 
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short distance region. The scale A plays the role of an UV cutoff. The logarithmic terms can 
be summed to all orders in QCD in the leading log approximation by applying RG-improved 
operator product expansion to (|1|) JTI], [L^, |13|]. Applying OPE to J^(x)J^(0), 

J J*{x)J*{U)D^{ X , A)dx = ]T Q(/i/A, a(/i))Of (0), (17) 



with q satisfying 



where 



fi 2 -^ + Piac)-^ + 7*0)1 Ci(ji/A, a(/i)) = 0, (18) 



7 j = 2 7j - 7j , (19) 
and 7j, 7, are the anomalous dimensions of the current and the operator Oi respectively, and 

d 

(3(a) = fx 2 —a = -p a 2 (l + O(a)), 
Po = ^(11 -jjivy), (20) 

fi(a) =7,«« + 0(a 2 ). (21) 

Note that 7^ arises only from the UV divergences in diagram (c)-(f). 
Integrating the RG equation, 

Cj(/i/A, a(/i)) = Q(a(A))exp / da , 



where 



«»0») /3(a) 

cf^Ca^.aCA)) (22) 



L 4 (aOi),a(A))=(^ "° . (23) 

The constants cj ' 1 can be determined from the tree-level amplitudes. Substituting (|2^) into 
(O), the one- loop effective Lagrangian for (JI|) should read 

C int (0)=J29 2 oV(»,Vc%0* (24) 



In computing the matrix element < f\Of-\i>, we take the MS subtraction scheme for the UV 
divergences. The scheme dependence in the matrix element is compensated by that of q to 
make physical amplitudes scheme independent. 

Since the EHLQ effective action at one-loop must assume the same form as in the tree-level 
action, the short distance effect in diagram (c)-(f) results in as a mixing among 77's with an 
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appropriate scaling by Lj. Let us first consider u^d^ —> u^d^ process. The effective Lagrangian 
for this process is 

£ UL d L = (2r/ - —)Oi + r?i0 2 

= (2r;o + |)0 + + (^-2^o)0-, (25) 

where 

Ox = u L Yu L d L ^d L , 

2 = u L ^d L d L ^u L , (26) 

and 

± = \{0 2 ±O l ). (27) 

Because + (0 ) is symmetric (anti-symmetric) under interchange between up and down 
quarks and separately up and down anti-quarks, 0± are multiplicatively renormalized to all 
orders in QCD, and their anomalous dimensions at one- loop are given by |TB 



^=^1^^ (28) 



According to (0), the effective action at one-loop is 

C ULdL - (2770 + f )£ + 0* + (^ - 2 Vo )L_0 R 



= \ ((2^ + |)L + - (^ - 2r ?0 )L_) O* 

+ i ( (2770 + ^)L + + - 2no)L.) O* (29) 



with 



L ± {a^),a{A))=(^r) * . (30) 



,«(A) 

Thus at one-loop order we have to replace 770, 771 by 770, 771 defined by 

( fjo\ ( \L + + \L_ UL+-LJ) \ ( 770 



(31) 



Let us now consider ^l^r — > mli^r which involves a color octet current. The Lagrangian for 
this process is 

£u L d R = VdOi + r] 8d 2 (32) 
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where Oj now are defined by 

O x = u L ^u L d R ^d R 

2 = u Ll »—u L d L Y—d R . (33) 

Unlike in the previous example, in this case there is no simple symmetry argument to find 
multiplicatively renormalized operators, and so we are going to diagonalize the one-loop mixing 
matrix explicitly. 

From the UV divergence part in (a)-(f), we have 

(£)'-*(£)"■ 

where 



(4tt) 2 e 

3(iVg-l)g 2 1 1 . 3(iV r 2 -2)g 2 l 
27V 2 (4tt) 2 e 1 /V c (4tt) 2 e 

with e = |(4 -n). The one-loop anomalous dimension of Oi is then 



(35) 



Diagonalizing T, we have 



with 



L = ( | 8 t ] . (38) 

V 27 9 / 



The one-loop anomalous dimension of octet current arising from (a) in Fig.l is 



Thus at one-loop level, 
where 



(!) Nc fon\ 

7 * = -to' m 



CuLdK^VdOF + fkiO? (40) 

(41) 



oil \ L 8 





(43) 



and 



^8± 



with 



Similar calculation gives 



and 





A 1 ! /p* 



-(7S-27«)/A> 



(44) 



J 1 ) 
X 1 ) 



-1 

47T 
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^ = (45) 



(46) 





(47) 



where L± are defined in (|30|). The transformation for 7? u , r/ 8u can be obtained by replacing r] d , r] 8d 
in ([41]) with r] u ,risu respectively. The modified coupling 77's should also be used in calculating 
the 2^3 tree-level amplitudes. For notational convenience, we keep using 77's instead of the 
modified couplings; however, in the rest of the paper, all 77's should be understood as 77's defined 
in m, (m, <m, and ©. 



5 Calculation framework 

When computing the one-jet inclusive cross section, we need the one-loop QCD amplitudes 
for 2^2, not the squared amplitudes, since QCD and the composite model interaction act 
coherently. To use the one-loop QCD helicity amplitudes for 2 — > 2 by Kunszt,et.al. || 
calculated in the 't Hooft Veltman scheme [14]], we also calculate the diagrams in Fig.l in the 't 



Hooft Veltman scheme. Kunszt,et.al. also isolated the soft and collinear divergences in 2 — > 2 
and 2—^3 processes, exposed explicitly the cancellation of these divergences among them, 
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and gave a complete prescription for the one-jet inclusive cross section @, 0. In this section 
we briefly review the calculation scheme of Kunszt,et.al., and identify terms to be calculated 
for the one-jet cross section. We follow the notation in || [7|] and readers should consult the 
references for more detailed discussions. 
The one-jet inclusive cross section 



d(Tjet 



(4£ 



can be written as 
where 



dpjdyj 

7(2 ->2) +7(2 ->3), (49) 



1(2-2) = ydp 2 da{2 d ^ 2) SM,pZ) 

/(2-3) = $JdP* M l^ 3) SM,&P£) (5°) 



and 



cfp 2 = dy 1 dp 2 dy 2 d(f)2 

dp 3 = dytdp 2 dy 2 d(p 2 dp 3 dy 3 d(f) 3 . (51) 

5*2, S3 define a jet algorithm, and pi, y.i denote transverse momenta and pseudo-rapidities of the 
partons respectively. 

7(2 — ► 2) can be divided into singular (~ l/e p ) and nonsingular parts 

7(2 -> 2) = 7(2 -> 2) s + 7(2 -> 2)^ (52) 

with 

7(2 - 2) NS = ^J dp 2 p 2 J2L AB (V (4 W) + ^S(a,p)) S 2 (pM, (53) 

where 

f A (a A ,x A )f B (a B ,x B ) (KA s 

l ab — \ 7 — s , (o4J 

w(a A )x A w{a B )x B 

and a = (a A , a B , Oi, ^2) for parton flavors, and p = (pa,Pb,Pi,P2)- Indices A, B and 1, 2 denote 
the initial state and the final state partons respectively, ijj^ is the Born amplitude squared 
and tpffs i s the nonsingular part of ip^ defined in 

J2 \A(a A + a B ^a 1 + a 2 )\ 2 =g 4 U 4 \ a ,p) + 2g 2 c r (-^-) V (6) (a, p) + 0(g 4 )) (55) 

colors \ X^ES/ I 

spins 
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where Qes is an arbitrary scale introduced by Ellis and Sexton 0J, and 

1 r(i-e) 2 r(i + e ) 



cr 



(4tt) 2 - £ T(l - 2e) 



(56) 



The singular part 1(2 — > 2)g depends only on the Altarelli-Parisi functions and the tree-level 
amplitudes in four dimensions, ^M 4 )(a, p) and ^£'(a, p), with the latter defined by 

& c) (a, p) = ^T^T^ J] ^41c lCs (2 - 2)4^0,(2 - 2). (57) 

For a = (qi,qj,qi,qj)i T a = X a /2 for the final state quarks and T a = — (A a )*/2 for the initial 
state quarks. 

Similarly 1(2 — ► 3) can be divided into singular and nonsingular parts 



J(2-3) 



J(2-3) 4(2 -3) +E4(2-3) 

n J n 

/ /inite (2^3)+^/;(2^3) 



(58) 



where n runs over A, 5, 1, 2. The soft and collinear divergences are contained in I' n (2 — > 3), and 
1^(2 — >■ 3) is divergent only when p 3 becomes soft or collinear to the parton n. J/j n j te (2 — > 3) 
is by construction well defined over all parton phase space, depends only on the tree-level 
amplitudes in four dimensions, and so the phase-space integration can be done numerically. 
Separating collinear divergence from soft divergence, I' n (2 —> 3) can be written as 



with 



/;(2^3) = C /i (2^3) + C«(2^3), 



C*(2-3) = I™ ft (2^3) s + I™ ft (2^3) NS 
C"(2-3) = C"(2 -> 3) s + C"(2^ 3)* s , 



(59) 



(60) 



where the explicit form of each term is given in |], 0. The singular and nonsingular terms in 
(|B0p involve Altarelli-Parisi functions and only tree-level amplitudes in four dimensions. For 
example, 



r 2 oft (2 -> 3) 



NS 



1ts o ; dp 2 p 2 J2 L AB [ij S 2 ° ft (^P)\ NS S2(Pl,P2] 



with 



ATS 



(4,c) ™ 
2m 1 2mi 



(61) 
(62) 



m=A,B,l 



where T 2m is a universal function of s^- = (pj +Pj) 2 - 
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Adding Q5%| ) into (f)2"D, we have 

/ = 7(2 -> 2) W5 + £ (C /4 (2 -> 3W + C"(2 - 3)jv S ) + J /infte (2 -> 3) (63) 

n 

with complete cancellation of the singular terms. From (|63|) we see that for the one-jet inclusive 
cross section, we need to calculate the tree-level amplitude V^j^mn^ ^ n f° ur dimensions, and 
the one-loop amplitudes ip^NS °f defined in (|55|). 



6 One-loop amplitudes 



The general form of the amplitude for qi{pi)qj{j>2) — > qk{p^)qi{Pi) to one-loop for each helicity 
channel can be written as 



where 



A{p\ + p 2 -> P3 + Pi) = 9 c (A 5 u 5 kj + B dijSki), 



A = A<®+g 2 A<» 
B = BW+cfB®, 



(64) 



and 



A(i) _ 4W 1 4W 
^1 — si-nnn T" S± C onti 



,(<) 
> cont • 



(65) 
(66) 



1 QCL> T" ^cont) -° — n QCD -°, 

i(») S(*) ot1 j J(0 S(i) 
^-QCDi n QCD dIlu ^con^-Ocor 

fermion contact interactions respectively, and c is a channel dependent spinor matrix element. 
The amplitude squared for each helicity channel is then 



qcdi^qcd aiiU ^-cont-, con are ^ ne tree an d one-loop amplitudes arising from QCD and the 



Y.\A 2 = /|c| 2 iV c 2 {(^ (0) ) 2 + ( s(0) ) 2 + ^A^B^ + 2g 2 [A^Re(A^ 

colors c 

Comparing (|67D with (|55|), we have ^^Vjvs in — > channel, 
= £ |c| 2 iV 6 2 [(^ (0) ) 2 +(5 (0) ) 2 + — 

spins ' ' 

E l« r 



(67) 



7/; (6) 



spins 



A^Re(A^) + B^Re(B^) + — (A^Re(B^) + B<°>ite(i (1) )) 



,(68) 
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where 



A« = 


( QeS 








A ^ , 


I c r 


NS 


m = 


( Qes 








[[ V 2 , 


1 c r 


NS 



(69) 



For the one-loop amplitudes A^J nt , B^lnt, we calculate the diagrams in Fig.l in the 't Hooft 
Veltman dimensional regularization scheme in which the spins and momenta of internal particles 
are defined in n dimensions, while those of external particles are defined in four dimensions. 
The calculation of the one-loop diagrams in the 't Hooft Veltman scheme is much simplified 
since we can treat the Dirac 7 matrices as if they were defined in four dimensions. To show 
this, let us consider, as an example, the diagram (a) in Fig.l for UL(pi)d L (p2) — > u^p^d^pz). 
The diagram (a) is proportional to 



d n k 



<4-| 7 ^+ J* 4 )7 a (l-75)(A+ #)7ll-><3-|7 a |2-> 



2 / dxx dy d n k 



k 2 {k + Pi) 2 (k + pi) 2 
<4-|7 M (^+ |f)7 a (l -ls)U+ ^)7 M |l-><3-|7 a |2-> 



(70) 



where p = —xypi + (1 — x + xy)p 4i , and q 
n-dimensions, as 

r = iU + y\ 



(k 2 - x 2 y(l - y)s u ) 3 
1 — xy)pi — x(l — y)p4 L . Writing 7^, defined in 

(71) 



where 7^ are the four dimensional Dirac matrices and 



r 



for fx < 4, 



[ 7 5,7l = 0, for /i>4, 
the numerator in the integrand in (|70| ) can be written as 

< 4-1^(^4)+ P)% } U+ ^(4))7f 4 )l 1 -><3-|7r4)|2-> 
+ <4-|r kifa il7ll-><3-|7( Q 4)|2-> 



(72) 



(73) 



The second term in (|73|) is 0(e 4 ) and so it can be safely discarded because the soft and collinear 
divergence is at most 0(l/e 2 ). The four-dimensional Dirac algebra then gives the divergent term 



to 



<4-|7?4) A4)7 ( Q 4) A4)7f4)|l-><3-|7 ( Q 4)|2- 
8k 2 



> 



n 



-[12] <34>, 



(74) 
(75) 
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where k 2 is defined in n-dimensions. Then the integration over k in n dimensions can be done in 
the standard way |H|. For other diagrams we can similarly check that only the four dimensional 
7-matrices contribute to the loop diagrams. 

As mentioned before, the counter terms (g),(h) are nonvanishing only for the color octet 
currents. For the octet currents, from the fermion self-energy diagrams and (a),(b), they are 
given by 

~Y g2 (e) Cr ' " 4 * ree ' ^ 
where Atree denotes the tree amplitude of the corresponding contact term. Also the UV diver- 
gences in (c)-(f) should be subtracted in the MS scheme. For each helicity channel, we give 
, c, AW , i« , , BW in appendix. 



7 Drell-Yan Process 

If composite quarks and leptons share common constituents, exchanges between quarks and 
leptons of their common constituents would give rise to quark-lepton contact interactions at 
low energies. Then the signal from these contact interactions may appear in Drell-Yan processes. 
In this section we write down a general effective quark-lepton contact interactions and consider 
their one-loop QCD corrections in Drell-Yan processes. 

As in the EHLQ lagrangian, we consider a single family of quarks and leptons. Including 
more fermion families should be straightforward. With the first generation of fermions, the 
most general, helicity preserving, SU(3)qcd x SU{2) l x U{1)y symmetric quark-lepton contact 
interactions are 

Cql = ^{^qLYqLhriL + ^qL^Y^ lL ^Y lL 

HuIlI^IlUr^ur + idLl^Ldta^dR + £ e qLl^qLe R ^e R 
+iueU R ^u R e R Ye R + idedRl^d R e R ^e R + (£ s q l L d R e R V L + /i.e.) } , (77) 

where II = (v,e)L, q^ = (u,d) L . The last term in (|77|) is due to scalar (or pseudo scalar) 
exchanges and all the other terms are due to vector (or axial vector) exchanges. 

For the massless fermions, no amplitudes with the same fermion helicities in the scalar 
exchange term arise in the standard model, and so the contact term of the scalar exchanges 
provides the leading amplitude in that helicity channel. Therefore we may keep the amplitudes 
in the scalar channel at tree-level, and consider one-loop QCD corrections only in the vector 
channels. 

To calculate the one-loop QCD corrections in the quark sector in Drell-Yan process, we 
must consider virtual corrections in qq' — > IV and the real gluon emission in qq' — > IV G along 
with qG — > q'W . Since these processes occur only in s-channel in the lepton momenta, the 
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amplitudes for a given helicity channel factorize into a flavor-independent part and a flavor- 
dependent propagator part that also includes the couplings on the quark and lepton vertices. 
The contact interactions thus modify only the propagator part, and so the one-loop QCD 
corrections in this model are essentially identical to those in the standard model. This allows 
us to write the cross section at parton level in terms of the corresponding cross section with a 
virtual photon exchange in the standard model 



da W (qq' -> W) = da { ^(qq -> 11) 



D ( - h \qq' -> ll')q< 



QqQl 



(78) 



where h, Q q> i denote the helicities and charges of quarks and leptons respectively and q 2 is the 
invariant mass squared of the leptons. The helicity independence of dafi /dQ 2 , where q 2 = Q 2 , 
of virtual photon exchange allows us to write ([78]) as 



dQ 2 

Similarly for the qG — > q'W 



da^ da * 



dQ 2 



(qG - q'W) 



dQ 2 



da^ 
~dQ'- 



D^\qq' -> ll')q 2 



QqQl 



(79) 



(qG^qll) 



D^(qq' -> ll')q 2 



QqQl 



(80) 



From dcr 7 * /dQ 2 in Altarelli, Ellis and Martinelli we finally have the Drell-Yan cross section 
for IV pair production, 



da 



DY 



1 1 



dQ 2 



1 dx-i f 1 dx 



x 2 ) + (1 



436nsQ 2 Jo x x Jo ^{['/'^ 

6(1 -z) + a s (Q 2 )9(l - z)(f q>DY (z) ~ VM) 



+ 



[(qf(x x ) + qf(x 1 ))G^(x 2 ) + (1 <-> 2)] a s (Q 2 )0(l - z)(/ G ,w(^) - / G , 2 (*))} 



5i) 



with 

/G,fly — /g,2 



2 

3tt 
1 

47T 



1-2 



6-4z + 2(l + z 



2 . , ln(l - z) 



1-z 



- - 5z +-z 2 + (z 2 
2 2 v 



;i-^) 2 )ln(l-z) 



[1 + 1^)5(1- z) 



i2) 



where z = Q 2 /x\X2S and a/s is the invariant mass of the incoming hadron system. 
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From the standard model lagrangian and ([77|) we can read off the nonvanishing propagator 
part 

D(u L d L 
D(u L u L 
D{u L u L 
D(u R u R 

D(u R U R 

D{d L d L 
D{d L d L 
D{d R d R 
D{d R d R 

where 

m_ = 2 — | sin 2 9 W ; u + = — |sin 2 ^, 
d- = — 2 + | sin 2 9 W ; d + = Isin 2 ^, 
Z_ = -2 + 4sin 2 ^; /+ = 4sm 2 8 w , 

and g is the SU{2) L coupling constant. 

Finally we would like to add a comment on calculating da DY / dpidyi in the framework of 
Kunszt, et.al., where pi,yi are the transverse momentum and rapidity of a designated lepton 
respectively. Since the essential difference of Drell-Yan process from the one-jet inclusive pro- 
duction is that the soft and collinear divergences in Drell-Yan process occur only in the initial 
state, all the formulae for the inclusive jet cross section can be used with only minor modifica- 
tion. First, since the final state leptons in Drell-Yan process are identifiable, the jet algorithm 
in (|50|) should be replaced by 

S 2 = 2\5{ Vl - Pl )% - yi ), S 3 = 3!<J(p, - - yi ), (84) 

and since soft and collinear divergences occur only in the initial state, the indices m,n in (p7|), 
(|58|), and (|63|) should run only over the initial state. Also the restriction on p 3 that it is the 
smallest among the final state parton momenta should be revoked. 



v L e L ) 2(?2 _ M ^ + 4A2 , 

— \ / 9 \ 2 u-l- QuQe . gfeo gfei 
eieL) = 1 - — ) — — H 5 h tttt _ tttt> 



Acos6 w J q 2 -M 2 z q 2 2A 2 8A 2 
9 V uJ + Q u Q e gli e 



' l!< " ( 1 cos 9 J q 2 - Ml + q 2 + 2 A 2 ' 

- , , 9 \ 2 u + l - . QuQe . gUu 



Acos6 w J q 2 -M 2 q 2 2 A 2 ' 

g \ 2 U + l + QuQe g 2 £,ue 
ene R ) = - — -= —r H =— + 



4cos6 w J q 2 -M 2 q 2 2 A 2 ' 

■> ! 9 \ 2 d_l_ Q d Q e g 2 £ g 2 ^ 

' L< ' ' " 1 4cos#J q 2 - M 2 q 2 2A 2 8A 2 : 

g \ 2 dJ + Q d Q e g 2 ^ e 
e R e n) = T—Tn 75. 777 + — 5- + 



Acos6 w J q 2 -M 2 z q 2 2A 2 ' 

g \ 2 d+L Q d Q e gli d 
e L e L ) = 7T -5 777 + — 5- + 



4 cos ^7 q 2 -Ml q 2 2A 2 ' 

/ 9 \ 2 d+l+ , QdQe . go£de /qq\ 

eR6R) = [l^KJ c^W 2 + + W (83) 
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For the virtual corrections, we need to consider only the diagram (a) in Fig.l, and its 
calculation goes exactly as in the case of the one-jet inclusive cross section. Details of the 
calculation and numerical analysis of the Drell-Yan cross section along with the one-jet inclusive 



cross section will be published elsewhere [16 
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A Appendix 



In this appendix we give the tree amplitudes c i ^conti B^Jnu an d the one- loop amplitudes 

A^ nt B^ nt defined in (0), (|64]),(|65|) and (^J), and for completeness, corresponding terms in 
QCD for the following helicity channels: (i) u^d^ — > ulcLl, (ii) ulUl — > ulUl, (iii) u R d^ — > 
u R d L , (iv) u L d R -> u L d R , (v) u R u R -> wrWr, (vi) -> and (vii) u fl dji -> M fl d fl . 

All the momenta of external fermions are assumed to have positive energies, and Sy = (p?.+Pj) 2 - 
c : 

c ici 



(Al) 



A^.B^: For contact terms, 



in, 
(iii 

(iv 
(v 



12 



(i) 


-2i[12] <34> 


4 s 2 


(ii) 


-2i[12] <34> 


4 s 2 




-2i[24] <13> 


4 s 2 


» 


-2i[13] <24> 


4s 2 


(v) 


-2i[34] <12> 


4s 2 


[vi) 


-2z[34] <12> 


4s 2 


vii) 


— 2t[34] <12> 


4s 2 



4(0) 0(0) 

-^cont ^cont 



i) — 2r?o + ?7i/2 -77! 

-2(r?o + V4) -2(770 + 771/4) 

-?7« + V8u/2N C -r)8u/2 

-rjd + Vsd/2N C -r]8d/2 

(vi) -2r/ dd -277^ 

(vii) ~Vud -vL 

and for QCD, A^ CD Bf CD 

(i) -l/2N cSli l/2si4 

(ii) -l/2iV c s 14 + l/2s 13 -l/2iV c Si 3 + l/2si4 

(iii) -l/2iV>i4 l/2s 14 

(iv) -l/2N c s u l/2s M 

(v) -l/2iV c 5i4 + 1/2S13 -l/2iV c 5i3 + 1/2S14 

(vi) -l/2iV>i4 + l/2si3 -l/2JV c s 13 + l/2si4 

(vii) -l/2iV c si4 l/2si4 



(A2) 
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tl>mn'- For each helicity channel, 



^(4) 

4i c) 



( A (o)) 2 + ( B (o)) 2 + ^(o )jB (o) 

(AW) 2 +(5(°)) 2 + 2N C A^B^ 



\~c\ 2 N* 

^/N c - \c\ 2 N c 
-ip ( ^/N c +\d\ 2 N, 
-ipW/N c +\c\ 2 N, 



4r\ 4r 



A (o) 



2 / ,„>\2 



where #),5 (0) 

are the sum of the corresponding contact and QCD terms as defined 
A (1) , B (1) : With the auxiliary functions, 



F 3 



2(N 2 - l)/N c 



3 + fin 



S14 



• • I • ^ 111" 



•514 



3 
2 

F 2 (.Sl2 



4 + fin 



^ES 



and 



lo 



2o 



-In 2 

l/N e 
+ ln 2 



S13 

-sis) 



1)/N 2 



| (In 2 



4+ fin 



S13 



In 2 



+ |ln 



J 2- 



3es 
-bis 



i In 2 



— S12 



3gg 

SIS 



5-iV? 



+ |ln 



| In 



^ + 3 In 




the contact terms are 
(i) 

(ii) 
(iii) 

(iv) 

(v) 
(vi) 



.4, 

A 
B 
A 

B, 

Replace 
Replace 
Replace 
(vii) Replace 



(i) 

cont 

(1) 
cont 

(1) 
cont 

(1) 
cont 

(1) 
cont 

(1) 
cont 



l-m fin 2 



(2r/ - (F x - F 2 /iV c ) + ^(pg <-> p 4 ) 
(2t/ - ^/2)F 2 + 77! (Fx - F 2 /N c )(p 3 <-> p 4 ), 
2(770 + 77i/4) [F - F 2 /N c + F 2 (p 3 <-> Pi)] 
2(770 + 77i/4) [F 2 + (Fx - F 2 /N c )(p 3 <-> p 4 )] , 
77«(Fi - F 3 /iV c ) + 77 8u (F lo - F 2o /iV c ) 
77«F 3 + 77 8n F 2o , 



Vu -> ?7rf, ??8u 
(7/0 + ^/4) - 
(7/0 + 7/1/4) - 
(2770 - W2) " 



-> 7/ M m (in), 
?7 UU in (ii), 
Vdd in (ii), 

-> Vud, 771 -> 7/; d in (i). 



For QCD, from the one-loop amplitudes in 

(i),(vii); 

2(i) _ _ 

^-QCD — 

5(1) 
n QCD 



1 



2514^, 



#1 



lHl(^l + ^) 
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1(1) 

A QCD 

5(1) 
n QCD 



S14 

1 
Sl4 



//l/iV c +^[(^l + ^2)(p3^P4)] 



1(1) 

6(1) 
n QCD 



1 



2si 4 JV« 



2^(^1 + ^), 



iV c (13/9 + 7r 2 )-fiV / + 



S14 



+ ^(1 + Sl3 / Sl2 ) (in 2 



-Sl4 
S12 , 



-3N C + ^iV c 



2 In 



In 2 



(iV c -l/iV c )+47rA)ln 



+si4/si 2 ln(f^)] +ln 
+3/iV c + 2/iV c ln(^) 

-(iV 2 - l)/N c [1^(1 + Sl3 / Sl2 ) (in 2 (g) + vr 2 ) 
W*i 2 lnfe)+21n(tf) Infe)" 



+ 7T 2 ) 



-^14 
S12 



3gS 
S14 



7V c (13/9 + vr 2 ) - f ^ + 8/N c + (N c + 1/N C ) ^(1 + s 12 /s 13 ) 



(ln 2 (^) +7 r 2 )+ Sl4 M3ln(^f) 



In 



-2N r In 



+ ln 2 



-Si4 
S12 



- |A^/ + 3/A^ c + 2/iV c ln 



Q 



ES. 



S14 



;-iV c + l/iV c )+47r/3oln 



S14 
-«13 



-3iV c + ±±iV c 



ES 



-(iV 2 - l)/iV c [^(1 + Sl2 / Sl3 ) (in 2 + vr 2 ) 

W*i3ln(^)+21n(%)ln(^) 
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Figure Captions: 



Fig. 1: One-loop diagrams from contact interactions, (g), (h) are the counterterms for the 
octet current renormalization. 

Fig. 2: One-loop diagram from dimension six operator, (b) is suppressed by a factor s/A 2 
compared to the diagrams in Fig. 1. 

Fig. 3: Penguin diagrams induce form factors in quark-gluon vertex. 
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(e) 






(h) 



(f) 



(g) 



Fig.l 




Fig.3 



